Collective Behavior of Asperities in Dry Friction at Small Velocities 
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We investigate a simple model of dry friction based on 
extremal dynamics of asperities. At small velocities, corre- 
lations develop between the asperities, whose range becomes 
infinite in the limit of infinitely slow driving, where the system 
is self-organized critical. This collective phenomenon leads to 
effective aging of the asperities and results in velocity depen- 
dence of the friction force in the form F ~ 1 — exp( — l/ii). 
PACS number(s): 46.30.Pa; 64.60.Lx 



I. INTRODUCTION 

Phenomena connected with mechanical properties of 
complex systems are subject of intensive study in the last 
decade. Generally speaking, the difficuhy stems from the 
fact that both macroscopic scale and mesoscopic scale are 
important. For example, the contact area of two grains 
of sand is a mesoscopic object, but its properties result in 
macroscopic behavior of a sand heap. Among the whole 
family of such problems, the dry friction emerged in re- 
cent years as a hot subject. Besides the intrinsic interest 
in dynamics of contact interfaces sliding on top of each 
other, there are various systems studied recently, in which 
friction forces are dominant interactions determining the 
behavior. As examples, we may note two notoriously 
known phenomena: sand heaps and earthquakes. Equi- 
librium stress distribution in heaps of granular materials 
exhibits complicated localized structures . Dynamics 
of tectonic plates gives rise to power-law distribution of 
earthquakes, formulated in Gutenberg- Richtcr law [^,1). 
A one-dimensional counterpart of friction is e. g. the 
dislocation movement, which is responsible for the plas- 
ticity of metals. 

At least three regimes of friction may be distinguished. 
First, dry friction corresponds to tangential force acting 
on the contact of two macroscopic solid bodies. The slot 
between the bodies is empty. The friction emerges as a 
result of the rheological properties of the sliding bodies 
both at the macroscopic and mesoscopic scale. Second 
case, the lubricated friction, differs in the fact, that the 
slot between bodies is filled with a liquid and the mechan- 
ical properties of mesoscopic portions of the lubricant are 
responsible for friction. Third, friction of a single micro- 
scopic tip on a surface may be measured, which explores 
the microscopic properties of the surface Here we 
concentrate on the first possibility, dry friction. 

Dry friction is intensively studied today, both experi- 
mentally and theoretically P,|7|] . The commonly accepted 



picture of friction is based on the dynamics of a system 
of mesoscopic contacts, called asperities, scattered on the 
surface of the sliding bodies ^-p^ . The typical size of as- 
perities is constant, while their number is proportional to 
normal load. Hence the Amontons-Goulomb law, stating 
that the apparent contact area of two macroscopic bodies 
does not matter. 

However, many features are not well understood, e. 
g. the velocity dependence of the friction force. It 
is explained either as a consequence of the plasticity of 
the asperities, which is considered as a thermally acti- 
vated process (this phenomenon is called aging of the 
asperities) or purely geometrically, based on self-afhne 
shape of the surfaces Within the approach based 

on plasticity, logarithmic dependence of the age of the 
asperity on time is supposed on the basis of experimen- 
tal data, which suggest logarithmic velocity dependence 
of the friction force. On the other hand, geometrical ap- 
proach gives friction force proportional to v^^ for large 
velocities, while the behavior for small velocities depends 
on the fractal geometry of the surface. 

In the description of the process of friction two levels 
may be distinguished. On the global level, the averaged 
effect of asperities can be successfully described using the 
elasto-plastic model developed by Garoli, Nozieres and 
Velicky [||j9). This approach is effectively a single-site 
one. Only one asperity is changing its state and effect of 
all other asperities is described by effective surrounding 
medium. The scheme resembles remotely the coherent 
potential approximation (GPA) for electronic structure 
of alloys. The spatio-temporal correlations are consid- 
ered to be of very short range, and the mutually sliding 
surfaces behave in uniform way. 

On the other hand, the local, short time level of de- 
scription must take into account processes which hap- 
pen at several (or many) asperities simultaneously, or 
within a very short period of time, so that they can- 
not be considered as uncorrelated. Several approaches in 
this direction were already proposed, based on geomet- 
rical considerations pO| , pl| , on Frenkel-Kontorova [ p^ , 
Burridge-Knopoff and train models [ p"4HT^ or on an ex- 
tremal dynamics model with elastic interactions p^ . 

The extremal dynamics (ED) models are very appeal- 
ing, because they may grasp the "skeleton" of the prob- 
lem, despite their simplicity and rudimentary nature. 
Generally, ED is based on the assumption that only one 
site is evolving during one time step, namely the site 
which has maximum (or minimum: it depends on the 



1 



model in question) of the dynamical variable determin- 
ing the state of the system. However, the price to pay 
is that the time scale fixed by the frequency of the up- 
dates of single sites is not directly related to the real time 
measured in an experiment. 

Extremal dynamics models were successfully used in 
modeling various systems, like invasion percolation, bi- 
ological evolution earthquakes or dislocation 
movement [^,^. The model we propose here is based 
on the ideas of ED models, adapted to the fact, that in 
friction we are interested in macroscopic movement with 
non-zero velocity, while most of ED models are appropri- 
ate to the case of infinitely slow movement. 

Shortly, he evolution of our model proceeds at the most 
susceptible asperity, namely the asperity which bears 
maximum stress. A small mechanical perturbation, like 
release of stress at single site, may result in a burst of 
activity of large spatio-temporal extent. Following the 
terminology used in the theories of self-organized critical- 
ity , we will call such spatio-temporal areas of activity 
avalanches. The correlations present in the model will be 
described through statistical properties of the avalanches. 

From time to time, the ED of asperities is interrupted 
by a macroscopic "slip" of the body as a whole, in which 
all asperities are completely renewed. By combination of 
ED evolution with such macroscopic slips, we introduce 
non-zero macroscopic sliding velocity into the model. 

The rest of the article is organized as follows. In the 
section II the model is defined and the interpretation of 
the model parameters is given, in sect. Ill the presence 
of self-organized criticality is investigated for in the case 
of zero macroscopic velocity, while the effect of non-zero 
velocity on the breakdown of SOC as well as the veloc- 
ity dependence of the friction force is investigated in sect. 
IV. The section V summarizes the results and draws con- 
clusions from them. 



II. THE MODEL 

We propose the following model. There are N point 
contacts, asperities, each with stress b. The quantity b 
will be interpreted as the elastic energy stored in the as- 
perity. The model is one-dimensional (the generalization 
to realistic two-dimensional case is straightforward) with 
periodic boundary conditions , so the points form a closed 
ring. In each step, the point with highest stress bmax 
is found and released. The release of the stress means, 
that the point is removed. However, in order to keep 
the number of points constant, new point is introduced 
somewhere in the system. 

As a zeroth approximation, the location of new point 
may be chosen at random. However, in reality the po- 
sition of the new contact is determined by the detailed 
structure of the the surfaces of slide and track. The new 
contact is established at such place, where the surfaces 
are closest one to another. So, another number, d is at- 



tributed to each point representing the width of the slot 
between the surfaces, waiting in the vicinity of the as- 
perity for further updates (the actual slot directly on the 
asperity is zero, of course), and in the update the location 
of minimum d is found. Here a new asperity is reintro- 
duced. The values of b and d of the neighbors of the old 
and new sites are also updated. Generally, each site has 
K — 1 neighbors which are affected. For simplicity, we 
assume K = 2, and update only one neighbor (the right 
one). 

Let us allow for very slow motion of the slider as a 
whole. The energy stored in the released asperity may 
be transferred entirely to other asperities, or a part of 
it may be converted into kinetic energy E. It may also 
happen that some of the kinetic energy is returned back 
to elastic energy of some asperities. It is natural to ex- 
pect, that at higher velocities, the number of asperities 
affected by the transfer of the kinetic energy to the elas- 
tic one will be larger. We simplify this dependence by 
saying, that for E < i?thr only the nearest neighbors are 
affected, while for E > E^thr the slider slips macroscopi- 
cally over average distance XsUp. The average duration of 
the slip is Tgiip and after that time all parameters b and 
d of all asperities are newly attributed at random and 
E is set to 0. Then, the dynamics starts again. In this 
process, the kinetic energy ~ -Ethr the system had before 
the slip is dissipated. This makes a difference with the 
theories of one asperity dynamics, where the energy is 
dissipated immediately after release of a single asperity. 
In our model we do not describe the processes which hap- 
pen during the slip, e. g. we do not examine the energy 
dissipated in course of the slip. Similarly, we do not cal- 
culate the physical velocity corresponding to the kinetic 
energy E during the ED evolution. So, we isolate only 
those contributions to the friction force and the macro- 
scopic slider velocity, which originate in the ED process 
interrupted by instantaneous slips. 

The average macroscopic velocity Av stemming from 
the slips depends on the average time interval between 
two subsequent slips. We may determine this quantity At 
in the time units of the extremal dynamics process. Its 
relation to physical time is not straightforward, but we 
suppose that this ambiguity affects only units, in which 
we measure time and not the general dependence of the 
friction force on velocity. Thus, we write simply 

Aw = l/At (1) 

which corresponds to taking the average slip length Xsiip 
as length unit and average time needed to update sin- 
gle asperity as a time unit. The contribution Av from 
the ED process is dominant if the time between slips is 
much larger than the duration of the slip. At ^ Tgiip (i. 
e. slips are instantaneous events) and the real length 
trvelled between slips during the ED dynamics xed is 
much shorter than the slip length, xed ^ a^sUp- 

The contribution Ai^tii to the friction force coming 
from this process is then proportional to the energy dis- 
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sipated in one slip. Because we are using arbitrary units, 
we identify 

AFfri = E'thr • (2) 

Let us now describe tfie extremal dynamics of the 
model more formally. The model consists of N sites con- 
nected in ring topology. Each site i £ {1,2,..., N} is con- 
nected to its right neighbor r{i). The state of the model 
is described by the set (E, bi, 62, ^at, di, c?2, cIn) and 
the function r{i) which describe the connectivity of the 
sites. At the beginning, E = and both bi and di are 
uniformly distributed in the interval (0,1). The updating 
steps are the following. 

(i) Find maximum stress fomax = niaxi(6i) located at 
site imax- Remember its old right neighbor ioid — ''(*max)- 

(ii) Find minimum slot dmin at site imin- 

(iii) Change of connectivity: The site z^ax is removed 
by cutting its links to left and right nearest neighbors 
and re-inserted between iuiin and the site next to it on 
the ring. It will have new right neighbor i,icw = ^(*max) — 
r(jmin), and then set r(iini„) = «max- 

(iv) Kinetic effects: Set 

E' = E + S ■ 6,„ax, &max = (1 ^ '5)6niax, 

Ai = (bu - h^je{bM - b^^J, A2 = {bu ~ h^^„)e{bM - 
If £" > Ai + A2, we set 

E = E'-Ai- A2, bl^^ = h^,, + Ai, = b,_ + A2. 
if not, we set 

E = 0, bl^^ = + E'/2, b[_ = h_ + E'/2. 

(v) Stress redistribution: For ri,r2 random numbers 
distributed uniformly in the triangle 0<rl<r2<lwe 
set 

^»„a. = '^l^max, Kli = &Ld+(^2-?-l)&max. ^»„o„ = ^Lw + 
(1 - ''2)&;„ax- 

(vi) New values of slots d are attributed to old and new 
neighbors as well as to site imax, taking random numbers 
uniformly distributed in the interval (0,1). 

(vii) If _E > i?thr, slip occurs, which means that E is set 
to and bi and di distributed uniformly in the interval 
(0,1). 

The kinetic effects and the slip involve several param- 
eters. First, the parameter 5 describes how much of the 
elastic energy tends to be converted into the kinetic en- 
ergy. This parameter can be interpreted as a measure of 
the spatial density of asperities: in the elementary pro- 
cess of freeing a particular asperity, more elastic energy 
may be converted into the kinetic one, if the density of 
asperities is lower. If (5 = 0, the kinetic effects are turned 
off. 

The parameter 6m is the limit, up to which an asperity 
can absorb a portion of kinetic energy and convert it back 
to elastic energy. It should be the property of the surface 
itself, without any resort to the load and velocity of the 
slider. If 5 = 0, the parameter 6m does not enter the 
model. 



The slip is determined by the parameter E'thr- In niore 
realistic description, it would be necessary to introduce 
the function R{E) which would count the number of sites, 
including the extremal site imax, which are to be updated, 
provided the kinetic energy has the value E. Here we take 
the simplest form R{E) 3 + (iV - i)e{E - E'thr)- Even 
this parameter should be the property of the surface, ir- 
respectively of the load and velocity. 

Finally, we comment on the interpretation of the quan- 
tity iV, the average number of asperities. We suppose 
that it may serve as a measure of the external load. While 
N does not depend on the apparent contact area of the 
slider and the track, the asperity density, which is hidden 
in the parameter 5, does. 

So, in order to conform with the Amontons-Coulomb 
law, we expect that the surface properties will enter the 
velocity dependence of the friction force through the pa- 
rameter 6m and combinations Ef^^/N and 5/N . We will 
see later that it is exactly the case. 
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FIG. 1. Distribution of stresses P(6) and maximum stresses 
Pmax(6max) for N = 1000, S = 0.01, 6m = 0.9. The en- 
ergy threshold is infinite ( full line for P{b) and dotted for 
f'max(6max)) and Ethv ~ 0.08 ( dashed line for P{b) and 
dash-dotted for Pmax(&max)). Number of steps is 10®. 

III. INFINITELY SLOW MOVEMENT REGIME 

Let us first investigate the case in which no slips are 
allowed, which can be expressed by limit value i?thr = 
oo. In this case, the macroscopic movement is infinitely 
slow. If the elastic energy could not be transformed into 
kinetic energy E, i. e. if S = 0, the model would be 
a slightly more complicated version of the Zaitsev model 
for dislocation movement [p^ , which is known to be self- 
organized critical. The criticality manifested by power- 
law distribution of avalanche sizes is due to infinitely slow 
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driving. It is natural to expect SOC also in our model 
for 6 = 0. However, even for S > the condition of 
infinitely slow driving, which means technically that only 
one asperity is updated at a time, is also satisfied and 
SOC is expected as well. 

We simulated systems of size N = 1000. The first 
quantity we measured was the probability distribution 
of the stresses, P{b) and maximum stresses -Pmax(femax)- 
The function P(6) is continuous up to a critical value 
b = be and then suddenly drops to zero, which is behavior 
common in SOC extremal dynamics models. The value 
of be depends on d. The typical behavior is shown in Fig. 
|l|for S = 0.01. 

A fingerprint of self-organized criticality is e. g. the 
scaling behavior of the forward A-avalanche sizes 



/'fwd(s) = .s-^9(.s|A-A,|i/-) 



(3) 



The A-avalanche starts when 6max exceeds the value A 
and ends when 5max drops below the value A again. The 
size s of the avalanche is number of update steps from the 
start to the end of the avalanche. For numerical reasons it 
is simpler to investigate scaling of integrated distribution, 
Pf^^(s) = ds Pfwd(s) from which the exponents r and 
a can be determined. 
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FIG. 3. Rescaled forward avalanche distribution for 
TV = 1000, S = 0.001, 6m = 0.9. The critical threshold is 
Ac = 0.519 and the scahng exponents r = 1.27 and l/cr = 2.6. 
Number of steps is 10*. The corresponding thresholds A are 
indicated next to the symbols in the legend. 
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FIG. 2. Rescaled forward avalanche distribution for 
N = 1000, 5 = 0. The critical threshold is Ac = 0.7475 
and the scaling exponents r = 1.28 and l/cr = 2.6. Number 
of steps is 10*. The corresponding thresholds A are indicated 
next to the symbols in the legend. 

The Fig. |2| and Fig. || show the data collapse which 
confirms the scaling of the form (||). Best collapse was 
obtained for the following values of the parameters: (a) 
for (5 = we have Ac = 0.7475 and r = 1.28 and l/cr = 
2.6. and (b) for 6 = 0.001 and bu — 0.9 we have Ac = 
0.519 and r = 1.27, l/cr = 2.6. 



There is a minor difference in the exponent t giving 
best fit for 5 = and 6 = 0.001. However, we believe 
that this difference is within numerical uncertainty of the 
results and the model belongs to the same universality 
class irrespectively of parameter S. 

By qualitative inspection of the quality of the data col- 
lapse for different choices of the exponents, we estimate 
the error bars. Thus, we finish with the following critical 
exponents of our model: 

T = 1.27 ±0.02, cr = 0.38 ± 0.02 . (4) 

The forward avalanche exponent t is greater than in 
the ID Zaitsev model [^H, but close to the Snep- 
pen interface growth modcl~ po| , |l6| . Another ID model 
to be compared is the CDW model of Olami pi| , and 
anisotropic interface depinning model of ref. |19[| which 
have however significantly larger exponent r. The clos- 
est universality class seems to be the one of the Sneppen 
model (r = 1.26), but the value of cr = 0.35 in this class 
is smaller than in our model. 

Whether this difference is due to finite-size effect or 
the two models are in different universality class cannot 
be stated with certainty from our present data. Instead, 
we would like to stress a structural similarity of the two 
models, which may explain the similarity of exponents. 
Contrary to usual interface growth models |2^ , the Snep- 
pen model is a non-local one. After a single growth event 
(deposition of a single particle), an unbounded sequence 
of further steps is performed in order to re-establish the 
single-step property of the interface. So, the range of 
interactions fluctuates during the evolution, according 
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to the actual configuration of the interface. Similarly, 
the Zaitsev model, like most of other extremal dynamics 
models is local in the sense that after finding an extremal 
site, its neighbors are also updated, while the range of 
neighborhood is fixed. On the contrary, our model, like 
the Sneppen model, does not have fixed range of interac- 
tions, but it is established by the position of the minimum 
of the quantity d (the slot). We simulated also a version, 
in which the site, where new asperity is inserted, is cho- 
sen at random, instead of using the slot d. In this case we 
observed mean-field behavior characterized by exponents 
r = 1.5, a = 0.5. 



IV. FRICTION AT NON-ZERO VELOCITY 

In the preceding section we dealt with stationary prop- 
erties of the model. In order to account for macroscopic 
movement, transient properties are of interest. First, we 
investigate the evolution of the kinetic energy E and its 
approach to the stationary value Eoo, if we forbid the 
slips, i. e. Ethr = oo. In Fig. ^we show the time evolu- 
tion oi E/N for different values of the model parameters 
6, 6m and number of asperities N. 



to N, which is in turn proportional to normal load. Thus, 
we recover the Amontons-Coulomb law for static friction. 

The approach of the kinetic energy to its stationary 
value is exponential, as is demonstrated in the Fig. 
This type of approach is directly reflected in the velocity 
dependence of the friction force, as we will see below. 
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FIG. 4. Time evolution of tlie kinetic energy per asper- 
ity. The parameters are as follows; A'' — 1000, 5 = 0.01, 
6m = 0.9 (full line) iV = 10*, 5 = 0.001, 6m = 0.9 (dot- 
ted line) N = 1000, S = 0.001, 6m = 0.5 (long dashed line) 
iV = 1000, S = 0.001, 6m = 0.9 (short dashed line). 

The most important observation is that the station- 
ary value Eoo/N depends on 5m, while the dependence 
on S and N is within the noise level. (We observe that 
both large N and small S suppress the relative fluctua- 
tions of the kinetic energy around the stationary value.) 
The physical significance is clear: the static friction force, 
which is according to eq. (^ equal to E^c , is proportional 



FIG. 5. Approach of the kinetic energy to its stationary 
value, for A*' = 10*, S — 0.001, 6m = 0.9. Stationary value is 
taken as Eoo/N = 0.155. 

If we set the threshold i?thr < Eoo, quasi-periodic be- 
havior is observed: the kinetic energy grows, until it 
reaches the value of the threshold, and then the system 
is reinitialized. This regime is illustrated in the Fig. 
If the threshold is close to E^o, the slips are less regular, 
due to fluctuations, but for smaller vales of the threshold, 
the slips occur with fixed frequency. The mean number 
of steps At between slips is determined by the way how E 
approaches the stationary value. Because Ethi is related 
to the friction force by (H) and mean period of slips to 
the velocity, according to (|l]) , the velocity dependence of 
the friction force is measurable in our model. The Fig. 
shows the results for various S and 6m- If we denote 
^0 = the static friction force, we observe by plotting 
the velocity dependence in semi-logarithmic scale that 
the following law is well satisfied 



AFm = Fo • (1 - exp(-A 



Ai;iV 



)) 



(5) 



with some constant A characteristic of the model. We 
have found A = 3.6 ±0.3. The deviations from the above 
dependence for Av < S/N are due to time fluctuations 
of E, which lead to less regular slips. However, as we al- 
ready mentioned, the relative fluctuations decrease with 
N, so we expect the dependence (^) to hold for all veloc- 
ities in thermodynamic limit N —^ oo. 



5 



0.14 - 
0.12 ■- 
0.1 - 
0.08 
0.06 - 
0.04 - 
0.02 





/I 







FIG. 6. Time dependence of the kinetic energy E, for 
= 10^, (5 = 0.1 and bM = 0.9. The slips occur in the 
moments when E drops to 0. 

As 6 is interpreted as a quantity proportional to the 
density of asperities, the ratio S/N does not depend on 

, i. e. on the normal load. Because also Fq was found 
to be proportional to N, the form of (|^) conforms with 
Amontons-Coulomb law for all Aw. 
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FIG. 7. Velocity dependence of the friction force, for 
N = 10^ S = 0.001, bM = 0.9 (+), 5 = 0.001, &a/ = 0.5 
{x),S = 0.01, bM = 0.9 (o). 

For large velocities the friction force decreases as 
Ai^fri 1/ Av. The same velocity dependence was found 
also using different approach |10|] . 

Now we turn to the influence of the macroscopic move- 



ment, connected to the slips on the self-organized critical 
behavior investigated in the last section. Each slip reini- 
tializes the values of b and d and the evolution towards 
the critical attractor begins from scratch. It means that 
the long-range correlations characteristic for the critical 
state cannot fully develop. The difference can be seen 
already in the distribution of stresses. Fig. I. The sharp 
edge in P{b) observed in the infinitely slow driving is 
smeared out. The position of the edge determines the 
critical threshold Ac for the forward avalanches, so we 
expect that no scaling of the type (^) will hold, as soon 
as the macroscopic movement has non-zero velocity. 
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FIG. 8. Distribution of jump lengths, for N = 10^, 
S — 0.001, bjvi = 0.9. The full line is case without slips 
(-Bthr = +oo). Dotted lines have slips allowed and the val- 
ues of Ethr/N are indicated next to the position where the 
lines reach the right edge of the figure. In the inset, distribu- 
tion of jump lengths, for S = and = 10^ (full line) and 

— 10^ (dotted line). Note that the inset makes it clear that 
the upward bent in the distribution for i?thr = -l-C)0 is mere 
finite-size effect. 

However, the most direct way how to investigate the 
breakdown of criticality due to the slips seems to us to 
be the calculation of the distribution of jump lengths. If 
in certain time step t the maximum stress was found at 
site it and in the next step at site it+i, we can compute 
spatial distance between this sites as follows. Let rt{i) be 
the function which determines the connectivity in time i, 
namely rt{i) is the site connected to i on the right-hand 
side. The jump length s is defined as follows: starting 
from if and applying rt we come to the right neighbor of 
the extremal site at time t, rt(it)- Then, applying s — I 
times the function rt+i we must end at it+i- So, s is such 
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that it+i = (rt+i'' o rt){it). 

In the self-organized critical state the probability dis- 
tribution of jump lengths is power law, Pjmp(s) = ■ 
For £^thr = oo it is actually observed in our model, as in- 
dicated in the inset in the Fig. |[ The comparison of the 
distributions for iV = 10^ and N — 10^ is shown in order 
to give idea of the magnitude of the finite-size corrections 
to the power-law behavior. 

The situation with non-zero macroscopic velocity, 
-Ethr < Eao, is shown in the main Fig. ||. When i?thr 
decreases, the velocity increases and the scale on which 
^jmp(s) obeys a power law shrinks. The correlations do 
not have time enough to develop on the scale of the whole 
system, but only on shorter distances. So, we may con- 
nect the velocity dependence of the friction force to the 
level of correlations between the asperities, which are 
present in the system. Contrary to the theories where 
the velocity dependence stems from the aging of a single 
asperity, here the aging is a collective effect. The age cor- 
responds to the range of correlations. For zero velocity 
the correlation length is infinite and the age is infinite as 
well. 



V. CONCLUSIONS 

We presented a model of dry friction based on the con- 
ception of slider and track interacting through a system 
of asperities. We proposed an extremal dynamics model 
in order to describe the processes during the movement 
of the slider. We found the decrease of the friction force 
with increasing velocity. For velocities approaching to 
zero, the friction force has finite limit. The origin of the 
velocity dependence is not in a change of properties of a 
single asperity, but in collective effects, involving many 
asperities. At zero velocity, the system is in a highly 
correlated, self-organized critical state. The value of the 
exponents are close to the Sneppen interface model, how- 
ever it is not clear from our data, whether the universality 
class is the same. 

Increasing velocity gradually destroys the correlations. 
It is possible to view the buildup of the correlations as a 
collective asperity aging mechanism, as a counterpart to 
the single asperity aging due to plastic deformation. Such 
collective aging leads to different velocity dependence of 
the friction force than in the models considernig single 
asperity aging and may be thus tested experimentally. 
Two-dimensional variant of our ID model would be nec- 
essary for real comparison. However, the generalization 
to arbitrary dimension is straightforward. 

This observation reveals also the limits of applicabil- 
ity of our model. It is appropriate in situations where 
the plastic deformation does not dominate. The model 
is applicable to the regime of very small velocities, where 
the usual logarithmic velocity dependence is inappropri- 
ate. It may be also used to describe friction over highly 
elastic surfaces, like rubber or some plastics, where the 



slow aging of single asperities may not be dominant. 
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